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摘 要 胡麻 产量 形成 模型 是 胡麻 生长 模型 中 的 重要 子 模型 ， 其 准确 性 关系 到 整个 模型 的 模拟 性 能 。 本 研究 
拟 构建 能 较 准 确 预测 不 同 生产 情况 下 胡麻 的 产量 形成 模型 ， 为 胡麻 实际 生产 提供 科学 指导 。 模 型 利用 
2012 一 2013 年 胡麻 西北 主 产 区 定西 和 榆 中 试验 站 不 同 肥料 、 不 同 播种 方式 、 不 同 种 植 密度 和 不 同 氮 磷 水 平 的 
试验 结果 ， 以 胡麻 生理 生态 过 程 为 主线 , 分 别 构建 了 产量 构成 因素 法 (基于 品种 遗传 参数 : 单位 面积 藉 果 数 、 每 
果 粒 数 、 粒 重 与 水 肥 胁 人 迫 因 子 、 累 积 光合 速率 ) 的 产量 形成 模型 和 基于 粒 壳 比 法 (基于 粒 过 比 和 薄 果 干 物质 总 量 ) 
的 产量 形成 模型 ; 利用 2014 一 2015 年 试验 区 实测 数据 对 模型 进行 较 广 泛 验证 。 验 证 结果 表明 , 产量 构成 因素 
法 对 定西 地 区 不 同 肥料 、 不 同 播种 方式 、 不 同 种 植 密度 产量 模拟 值 的 RMSE 平均 为 133.47 kg-hm“， 粒 壳 比 法 
的 RMSE 平均 为 195.51 kg:hm“; 利用 最 小 二 乘法 计算 的 决定 系数 及 产量 构成 因素 法 与 粒 过 比 法 分 别 为 0.819 8 
和 0.743 9。 榆 中 地 区 12 个 氮 磷 处 理 水 平 产量 构成 因素 法 平均 RMSE 值 为 140.70 kg.hm ”， 粒 壳 比 法 的 RMSE 
F 均 为 193.22 kg.hm-?; 产量 构成 因素 法 与 粒 壳 比 法 预测 的 产量 值 与 观测 值 R? 分 别 为 0.8329 和 0.805 8。 结 
表明 , 产量 构成 因素 法 对 产量 的 模拟 效果 优 于 粒 壳 比 法 ， 针 对 不 同 品种 不 同 栽培 方式 ， 把 单位 面积 薄 果 数 、 每 
果 粒 数 与 粒 重 作为 品种 遗传 参数 的 产量 构成 因素 法 能 更 精确 模拟 胡麻 的 产量 ， 具 有 较 高 的 预测 性 和 通用 性 。 
关键 词 ”胡麻 产量 形成 模型 ”产量 构成 因素 法 “ 粒 壳 比 法 “遗传 参数 
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Abstract The oilseed flax yield formation model was developed to predict biomass yield and provide scientific guidance for 
oilseed flax production, which was an important sub model of oilseed flax development model, APSIM-Oilseed flax. Its’ 
accuracy was related to the performance of APSIM-Oilseed flax. In the study, the oilseed flax yield formation model was 
developed for Northwest China using data generated from experiments, and was calibrated and validated under varying 
fertilizers, seeding methods, planting densities, nitrogen and phosphorus application conditions. The experiments were 
conducted at Dingxi and Yuzhong of Gansu Province from 2012 to 2015. Calibration was done using data for 2012-2013 and 
validation with data for 2014-2015. Two methods, yield component factors method and grain pod ratio method, were used to 
build oilseed flax yield formation model. The oilseed flax yield formation model of yield component factors was based on 
genetic parameters, such as, number of pods per unit area, grains number per pod, grain weight and water and fertilizer stress 


factors. The oilseed flax yield formation model of the grain pod ratio was based on grain pod ratio and total dry matter. The 
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results showed that the average RMSE for the model-simulated and field values related to the yield component factors method 
under different fertilizers, sowing densities in Dingxi was 133.47 kg per hm2. Then it was 195.51 kg per hm’ for the grain pod 
ratio method. Then R2 values were 0.819 8 and 0.743 9 for the yield component factors and grain pod ratio, respectively. The 
average RMSE was 140.70 kg per hm’ for the yield component factors method under varying nitrogen and phosphorus levels 
in Yuzhong and it was 193.22 kg per hm’ for the grain pod ratio method. The R* values were 0.832 9 and 0.805 8 of yield 
component factors and grain pod ratio, respectively. The validation results showed that the yield component factors method 
was better than the grain pod ratio method in simulating oilseed flax yield formation. In view of different varieties of oilseed 
flax under varying experimental treatments, the yield component factors method simulated oilseed flax yield with high 


accuracy and versatility using the number of pods per unit area, seeds number per pod and grain weight as genetic parameters. 


Keywords 


作物 产量 即 单位 土地 面积 上 作物 群体 的 产量 ， 
作物 产量 通常 分 为 生物 产量 和 经 济 产 量 , 一 般 产 量 
指 的 是 籽粒 收获 量 ， 即 经 济 产量 。 经 济 产量 占 生物 
产量 的 比例 ， 即 生物 产量 转化 为 经 济 产 量 的 效率 ， 
为 经 济 系数 或 收获 指数 。 通 常 ， 经 济 产量 与 生物 产 
量 呈正 比 ' 在 实际 生产 中 , 加 强 作物 生物 量 累积 动 
态 与 产量 形成 的 研究 有 利于 从 总 体 上 把 握 作物 生产 
过 程 中 。 因 此 ， 作 物产 量 的 形成 与 器 官 分 化 、 发 育 及 
光合 产物 的 分 配 和 累积 密切 相关 ， 了 解 其 形成 规律 
是 实现 作物 高 产 的 基础 趾 。 本 研究 则 在 通过 构建 胡 
麻 (Linum usitatissimum L.) 产 量 形成 模型 ， 以 期 对 胡 
麻 实际 生产 进行 科学 指导 。 

目前 ， 有 关 作 物产 量 形成 模型 的 研究 已 非常 广 
泛 ， 从 单一 性 模型 到 综合 性 的 模拟 模型 ， 从 一 种 作 
物 的 模型 到 多 作物 的 整合 模型 都 有 较 多 的 研究 ， 但 
已 有 的 产量 形成 模拟 方法 的 研究 存在 较 大 差异 。 目 
前 作物 产量 的 模拟 方法 主要 有 产量 构成 法 SM 收获 
指数 法 0 1、 分 期 转移 法 六 和、Logistic 方程 法 1 
以 及 全 球 气 候 变化 对 作物 产量 的 影响 模型 “和 采 
用 地 面 遥 感 和 高 空 融 感 观 测 数据 估算 作物 产量 等 模 
型 8°13。 邹 菩 等 中 在 大 麦 (Hordeum vulgare L.) 产 量 
构成 模型 研究 中 采用 产量 构成 法 构建 了 适用 于 不 同 
地 区 不 同 品种 的 大 麦 产量 模拟 模型 ; 郑 秀 琴 等 思 治 
用 小 麦 (7Tyiticum aestivum L.) 产 量 构成 因素 方法 ， 
建立 了 模拟 冬小麦 产量 形成 模拟 模型 ; 汤锅 等 所 
通过 粒 壳 比 和 角 果 物质 量 的 关系 计算 油菜 (Brassica 
campestris 工 .) 最 终 的 籽粒 产量 ; 张 亚 杰 忆 的 研究 中 
也 采用 了 衡量 角 果 本 身 干 物质 合理 分 配 和 库 源 关系 
的 重要 指标 一 一 粒 壳 比 计算 直播 油菜 的 籽粒 产量 ; 王 
新 等 3 构建 了 基于 收获 指数 的 加 工 番茄 (Lycopersicon 
esculentum Mill.) 产 量 预测 模型 ; 尹 红 征 等 中 和 郑 
国清 等 2 根据 玉米 (Zea mays L.) 籽 粒 干 物质 的 来 
源 构建 玉米 产量 形成 模型 ; 陈兵 兵 己 在 有 关 芭 麻 
[Boehmeria nivea (L.) Gaudich.] 产 量 模型 的 优化 研究 
中 ,运用 产量 构成 因素 法 、 多 元 线性 回归 、BP 神经 
网 络 和 GA-BP 神经 网 络 ( 采 用 遗传 算法 对 BP 神经 网 


Oilseed flax; Yield formation model; Yield component factors method; Grain pod ratio method; Genetic parameter 


络 进行 参数 优化 ) 分 别 建立 芝麻 产量 模型 ， 然 后 进行 
估 测 比较 ， 最 终 确 定 最 优化 模型 。 目 前 对 胡麻 产量 
模拟 模型 的 研究 除了 本 课题 组 在 文献 [26] 中 的 研究 
之 外 鲜 见 报道 。 胡 麻 产 量 形成 模型 是 胡麻 生长 模型 
中 的 重要 子 模型 ， 其 准确 性 关系 到 整个 模型 的 模拟 
性 能 。 胡 麻生 长 模拟 模型 是 依托 胡麻 的 生长 发 育 特 
征 在 计算 机 上 模拟 胡麻 生长 发 育 、 器 官 建成 和 产量 
品质 形成 等 与 环境 之 间 的 数学 关系 ,主要 包括 物候 
期 子 模型 、 叶 面积 指数 子 模型 5、 光合 生产 子 模 
型 (后 续 发 表 )、 器 官 建成 (后 续 发 表 ) 以 及 产量 形成 子 
模型 等 。 

本 研究 在 现 有 作物 生长 发 育 模 型 研究 的 基础 
上 ,以 胡麻 的 生理 生态 过 程 为 主线 , 分 别 采 用 产量 
构成 因素 法 和 粒 壳 比 法 建立 胡麻 产量 形成 模型 , 通 
过 估 测 比较 各 模型 的 模拟 精度 ， 最 终 确定 能 更 精确 
预测 胡麻 产量 的 模拟 方法 ,为 胡麻 的 实际 生产 提供 
引导 。 
1 材料 与 方法 
1.1 试验 设计 

试验 1: 肥料 试验 , 试验 因素 为 油 渣 、 磷 酸 二 
铵 、 复 合肥 。 试 验 设 不 施肥 作对 照 (CK)， 设 油 渣 施 
用 量 分 别 为 600 kg:hm“(Y1)、1 200 kg.hm ”(Y2)、 
2 400 kg.hm ”(Y3) 3 个 水 平 , 磷酸 二 铵 施用 量 分 别 
为 90 kg.hm (R1)、180 kg:hm™ (R2)、270 kg.hm 
(R3) 3 个 水 平复 合肥 施用 量 分 别 为 150 kg:hm” (F1)、 
300 kg:hm”(F2)、450 kg:hm“(F3) 3 个 水 平 , 共 10 
个 处 理 , 每 个 处 理 3 次 重复 , 共 30 个 小 区 。 小 区 面 
积 20 m”。 试 验 用 胡麻 品种 为 “ 定 亚 22 号 ’。 

试验 中: 覆盖 试验 ,， 设 3 个 试验 处 理 , 分 别 为 
残 膜 直 播 (T1); 残 膜 覆 至 春天 ,播种 前 揭 残 膜 ， 覆 
盖 新 膜 播种 (T2); 残 膜 覆 至 春天 ,播种 前 揭 残 膜 后 
直接 播种 (T3)。 各 小 区 氮 、 磷 、 钾 施肥 量 均 分 别 为 
112.5 kgN) -hm *\75 kg(P,0s)-hm * 52.5 kg(K2O)hm ; 
共 3 个 处 理 3 次 重复 。 小 区 面积 20 m”, 试验 用 胡麻 
品种 为 ' 陇 亚 10 号 *。 
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试验 川 : 播种 量 试验 , 设 7 个 播种 量 处 理 ， 分 别 为 
3x104 粒 ,hm ”(D1)4.5x10 粒 .hm (D2) 6x105 粒 -hm 
(D3)、7.5x104 粒 .hm (D4)、9x10 粒 .hm”? (D5)、 
1.05x10" 粒 .hm 了 (D6)、1.2x107 粒 .hm (D7); 各 处 理 的 
穴 距 和 行距 均 为 11 cm 和 15 cm, 每 穴播 种 籽粒 数 对 应 
种 植 密度 分 别 为 6 粒 、9 粒 、12 粒 、15 粒 、18 粒 、21 
粒 、24 粒 。 小 区 面积 13.3 m” 重复 3 次 。 各 小 区 气 、 磷 、 
钾 施 肥 量 均 分 别 为 112.5 kg(N)-hm*、75 kg(P2O3)hm-、 
52.5 kg(K2O)-hm ,胡麻 品种 选用 ' 陇 亚 10 号 *。 

试验 | ~ 川 于 2012 一 2015 年 在 甘肃 省 定西 市 西 
巩 驿 镇 (104"37'12"E, 3$"34'48"N) 进 行 。 试 验 区 属 时 
地 。 前 茬 为 全 膜 双 垄沟 玉米 ,前 一 年 玉米 收获 后 保 
护 地 膜 , 以 草木 灰 或 砂 土 覆盖 破损 处 , 冬季 避免 牧 
畜 践 踏 和 人 为 损坏 地 膜 ， 翌 年 春天 免 耕 直接 播种 。 
供 试 地 膜 为 聚 乙烯 吹 塑 农用 地 膜 ， 厚 度 0.008 mm， 
甘肃 省 天 水 天 宝 塑 业 有 限 责任 公司 生产 。 气 、 磷 、 
钾肥 品种 分 别 为 尿素 、 过 磷酸 钙 和 硫酸 钾 。 磷 、 钾 
肥 均 作为 基肥 施用 ; 氮肥 的 2/3 作为 基肥 ，1/3 作为 
追肥 于 现 蕾 前 追 施 。 小 区 间 走 道 30 cm, 重复 间 走 道 
50 cm， 四 周 设 1 m 的 保护 行 。 试 验 | 和 试验 中 的 播 
种 密度 为 750 万 株 .hm ,行距 15 cm， 穴 距 11 cm， 
每 穴播 10 粒 左右 种 子 。4 月 1 日 播种 , 8 月 5 日 收获 。 
栽培 管理 同 大 田 生产 。 

试验 IV: 氮 磷 施肥 试验 ,于 2012 一 2015 年 
在 甘肃 省 兰州 市 榆 中 县 育种 繁殖 场 (103"49'15"E~ 


104°34'40"E,35°34'’20”"N~36°26'20"N) 进 行 。 试 验 区 
属 水 地 。 试验 因素 为 所 和 磷 ， 采用 二 因素 随机 区 组 
设计 。 氮 设 3 个 水 平 ， 分 别 为 0 kg(N).hm (N0)、 
75 kg(N)-hm” (N11) 和 150 kg(N).hm ”IN2); 磷 设 4 个 水 
平 , 分 别 为 0 kg(P,Os)-hm” (P0)、75 kg(P;0s)-hm (P1)、 
150 kg(P20;)-hm“ (P2) 和 225 kg(P205)-hm“ (P3)。 共 
12 个 处 理 , 重复 3 次 , 共 36 个 小 区 。 各 小 区 均 施 钾 
肥 52.5 kg(K2O)hm-”。 所 、 磷 、 钾 肥 与 施用 方法 同 
上 。 各 小 区 灌溉 定额 均 为 2.7x10” m3.hm-“ (分 茎 期 灌 
1.2x10  m .hm “、 现 蕾 期 灌 1.5x10”m*.hm“)。 胡 麻 
品种 选用 ' 陇 亚 杂 1 号 :， 种 植 密度 为 900 万 株 .hm 
条 播 ， 播 深 3 cm, 行距 20 cm。 小 区 长 Sm、 宽 4 m， 
面积 20 m”。2012 年 3 月 24 日 播种 , 8 月 4 日 收获 ， 
2013 年 4 月 15 日 播种 , 8 月 9 日 收获 ,2014 年 4 月 1 
日 播种 , 8 月 5 日 收获 , 2015 年 3 月 25 日 播种 ,8 月 5 
日 收获 。 和 栽培 管理 同 大 田 生 产 。 

定西 地 区 地 处 黄河 中 游 黄 土 高 原 沟 到 区 ,海拔 
高 度 1 793 m， 年 均 气 温 7 'C, 年 日 照 时 数 2 500 h， 
无 霜 期 146 d, 年 降水 量 300~400 mm, 年 均 蒸发 量 
1 524.8 mm。 供 试 土壤 为 黑 坊 士 , 土壤 理化 性 状 见 表 1。 

榆 中 地 区 地 处 黄河 中 游 黄土 高 原 沟 饪 区 , 海拔 高 
度 1 793 m, 年 均 气 温 6.7 'C, 宇 10 ‘CC 积温 2 350 ‘Cd, 
年 日 照 时 数 2 563 h, 无 霜 期 146 d, 年 降水 量 300~ 
400 mm, 年 均 蒸 发 量 1 341 mm, 年 辐射 量 1 310 MJm。 
供 试 土壤 为 砂 壤土 , 土壤 理化 性 状 见 表 1。 


表 1 定西 和 榆 中 试验 站 土壤 理化 性 状 


Table 1 Basic physical-chemical properties of soil in Dingxi and Yuzhong experiment sites 


ee 合唱 yy ee 二 
试验 站 土壤 类 型 有 机 质 合 量 碱 解 氨 速效 砚 速效 钾 
Experiment site Soiltype Organic matter Total N Available N Available 了 Available K pH 
(gkg ) (gkg ) (mg'kg ) (mgkg ) (mgkg ) 
定西 Dingxi 黑 坊 士 Heilu soil 11.06 72% et 8.31 247.02 8.30 
榆 中 Yuzhong 砂 壤土 Sandy loam 16.56 59.01 13.83 127.62 7.75 


1.2 ”数据 来 源 与 测定 

模型 运行 所 需 参 数 为 气象 数据 (逐日 最 高 、 最 低 
气温 、 降 水 量 、 辐 射 ), 来 自 甘 肃 省 气象 局 农业 气象 
观测 站 ; 土壤 数据 (土壤 类 型 、 质 地 、 典 型 土壤 分 层 
饱和 含水 量 、 田 间 持 水 量 、 萎 六 系 数 、 容 重 、 有 机 
质 合 量 、 全 氮 、pH、 土 壤 碳 气 比 等 ) 来 自 农业 气象 观 
测 站 和 中 国土 壤 数 据 库 ; 品种 数据 (品种 类 型 、 发 育 
速率 、 干 物质 分 配 系数 、 比 叶 面 积 等 )、 田 间 管 理 ( 播 
种 期 、 播 种 深度 、 种 植 密度 、 施 肥 量 、 灌 溉 量 及 其 
时 间 等 ) 和 实际 产量 数据 来 自 试验 实际 观察 测定 。 

试验 观测 项 目 包括 胡麻 生育 时 期 , 各 生育 时 期 
各 器 官 根 、 芭 、 叶 、 果 的 干 物 质量 、 叶 面积 ,收获 
时 测定 产量 构成 要 素 及 产量 ， 播 种 前 和 收获 后 分 别 


测定 土壤 理化 性 状 及 土壤 水 分 动态 。 

1) 土 壤 基 础 参数 : 播种 前 分 别 测定 0~30 cm 土 
壤 的 铵 态 氮 含量 、 硝 态 氮 含量、 有机质 含 量 和 土壤 
合 水 率 ， 在 收获 后 测定 土壤 水 分 特征 曲线 和 容重 。 

2) 物 候 期 : 分 别 记 录 苗 期 、 现 蕾 期 、 盛 花期 、 子 
实 期 和 成 熟 期 。 

3) 叶 面积 和 干 物质 量 : 在 苗 期 、 现 蕾 期 、 盛 花期 、 
子 实 期 和 成 熟 期 进行 取样 和 测定 。 叶 面积 每 小 区 定 
株 10 株 进行 叶 面 积 的 连续 观测 ; 生物 量 每 小 区 随机 
取样 10 株 , 分 别 测定 荆 、 叶 和 戎 果 的 干 鲜 重 量 。 

4) 产 量 : 成 熟 期 按 小 区 进行 实 产 测定 ， 考 种 测 
定 产量 构成 三 要 素 。 

采用 SPSS 20.0 统计 分 析 软 件 进行 数据 整理 和 分 析 。 
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1.3 ”模型 调 参 方法 

本 模型 是 综合 胡麻 生长 发 育 模 型 (APSIM- 
Oilseed flax) 的 一 个 子 模型 。 模型 采用 VB.NET 设计 ， 
APSIM 已 建立 的 作物 模型 有 认 嘴 豆 (Cicer arietinum 
L.)、 绿 豆 [Vigna radiata (L.) Wilczek.]、 大 豆 [Glycine 
max (L.) Merr.]、 柱 花草 (Stylosanthes euianensias SW.)、 
花生 (4rachis hypogaea L.)、 人 蛋 豆 (Vicia faba 工 .)、 紫 
花 首 攻 (Medicago sativa L.)、 加 拿 大 油菜 (Brassica 
napus L.)、 小 麦 等 。 本 研究 通过 借鉴 APSIM-Canola 
加 拿 大 油菜 模型 0* 1， 根据 胡麻 生长 发 育 的 生理 生 
态 过 程 , 确定 影响 胡麻 产量 的 品种 遗传 参数 。 确 定 
作物 参数 过 程 : 首先 调查 收集 研究 区 作物 性 状 介绍 
与 试验 资料 ， 作 为 估计 作物 品种 资料 遗传 参数 的 基 
本 依据 ; 其 次 , 采用 基于 神经 网 络 的 投影 寻 踪 自 回 
归 BPPPAR(projection pursuit auto-regression based 
on error back propagation) 模 型 中 用 RAGA(real 
coded accelerating genetic algorithm) 优 化 投影 指标 函 
数 P ,调整 品种 遗传 特性 参数 。 
1.4 ”模型 检验 方法 

采用 模型 校准 与 验证 评价 中 常用 的 统计 指标 : 
均 方 根 误差 RMSE(root mean square error)、 决 定 系 数 
R*(determination coefficient) 和 相对 误差 REGelative 
error) 对 模型 进行 检验 ， 其 中 R? 用 于 评价 模型 的 预测 
能 力 , RE、RMSE 用 于 显示 模型 预测 中 的 误差 。 


[en 本 
RMSE = Za M5) (1) 
n 


> 0; — M)S;—S) 


2 _ 
| 本 Ci -MY NY 0 一 NO 


RE=|M-S/M (3) 
式 中 : S$, 和 M 分 别 为 模拟 值 和 实测 值 ，M 为 M; 的 平 
均值 , ”为 观察 值 的 数目 。 


2 ”模型 的 构建 

前 人 关于 作物 产量 预测 研究 有 多 种 方法 。 本 研究 通 
过 比较 产量 构成 因素 法 中 和 粒 壳 比 "两 种 方法 的 模 
拟 结果 ， 虽 在 找 出 能 更 精确 预测 胡麻 产量 的 模拟 方法 。 
2.1 产量 构成 因素 法 产量 形成 模型 

胡麻 产量 结构 由 单位 面积 上 的 浆果 数 、 每 果 粒 
数 和 粒 重 3 个 因素 构成 。 因 此 ， 本 模型 胡麻 潜在 产 
量 表示 为 : 


R (2) 


-6 
=N,xPxW,x10 (4) 
hh=BxD (5) 


式 中 : 及 是 胡麻 潜在 产量 (kg:hm“)，WN, 为 每 果 粒 数 ， 
PP 是 单位 面积 区 果 数 (个 .hm“)，W, 是 千粒重 (8)， 
只 是 单 株 闲 果 数 , DD 为 胡麻 群体 密度 ( 株 .hm )。 
单 株 作物 的 籽粒 数 (N,) 由 开花 期 茎 重 决 定 玉 ; 
N, = Re XW (6) 
式 中 : 到 是 开花 期 基干 重 (g)，R, 是 每 克基 的 籽粒 
数 (个 .g )， 取 值 133 个 .g 。 
因而 , 每 果 粒 数 表 示 为 : 
N,=N,/D, (7) 
胡麻 实际 产量 受 水 肥 胁 迫 因子 影响 ,水肥 万 
缺 因 子 取 水 分 和 氮肥 胁迫 因子 的 最 小 值 。 表 示 为 : 


y=, xmin(f,,f,) (8) 
/w= 工友 (9) 
CN 一 CN 
大 三 N Nmin (10) 
CNerit CNmin 


式 中 : 了 为 胡麻 实际 产量 (kg:hm“); ff, 、 是 水 分 
和 氮肥 亏 缺 因子 ; 双 为 胡麻 群体 实际 蒸腾 量 ; 为 
胡麻 群体 潜在 蒸腾 ,由 水 分 平衡 模型 计算 得 出 ;CN 是 
进入 叶 组 织 的 实际 氮 浓 度 ; CN sit、CNain 是 临界 氮 浓 
度 与 叶片 自由 生长 氮 浓 度 ， 由 氮 平衡 模型 计算 得 出 。 
土壤 水 分 平衡 、 氮 平衡 模型 由 本 研究 组 后 续 发 表 。 
2.2” 粒 壳 比 法 产量 形成 模型 

前 人 研究 表明 粒 壳 比 ( W/W,， 籽粒 重 和 果 壳 重 
的 比值 ) 是 衡量 靖 果 本 身 干 物质 合理 分 配 和 “ 库 源 ” 关 
系 的 重要 指标 之 一 “不同 品种 作物 粒 壳 比 有 显 
著 差 异 。 因 此 ， 本 研究 采用 粒 壳 比 计算 胡麻 的 籽粒 
产量 ,， 且 该 参数 为 品种 参数 ， 计 算 公式 如 下 : 


W, /W, 
y=C (11) 
时 1+W,/W, 


式 中 : 驴 为 籽粒 潜在 产量 (kg:hm“)，C, 为 胡麻 收获 
时 的 薄 果 干 物质 总 量 (kg-hm )。 
3 ”结果 与 分 析 
3.1 ”参数 确定 
胡麻 产量 形成 模型 参数 主要 包括 : 单位 面积 痕 
果 数 、 每 果 粒 数 、 粒 重 、 单 株 薄 果 数 、 氨 肥 胁 迫 因 
子 等 各 项 指标 , 如 表 2 所 示 。 
3.2 ”模型 验证 
3.2.1 ”产量 构成 因素 法 与 粒 壳 比 法 的 对 上 比 检 验 
分 别 采 用 产量 构成 因素 法 和 粒 壳 比 法 计算 胡麻 
产量 , 通过 选择 两 个 试验 站 相应 的 肥料 、 播 种 方式 、 
种 植 密度 和 和 气 磷 水 平 ， 调整 相应 参数 ， 运 行 模型 得 
出 产量 模拟 结果 ,其 RMSE 和 R? 值 见 表 3、 图 1。 
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表 2 胡麻 产量 形成 模型 品种 遗传 参数 


Table 2 Genetic parameters of oilseed flax yield formation model 


单位 面积 蒲 果 数 每 果 粒 数 粒 重 单 株 戎 果 数 单 株 籽粒 数 ”水 分 胁迫 因子 氮肥 胁迫 因子 ” 粒 壳 比 


品种 
Re Number of pods per Number of Grain weight Number of pods Number of grains Stress factor Stress factor of Grain/pod 
unit area (pods:hm *) grains per pod (g) per plant per plant of water nitrogen ratio 
i 
陇 亚 杂 1 号 8 148 7.60 7.50 15.52 118 0.75 0.68 1.05 
Longyaza 1 
陇 亚 10 号 11 655 8.20 8.13 22.20 182 0.75 0.68 1.12 
Longya 10 
| 可 
定 亚 2 10 920 8.15 7.43 20.80 170 0.75 0.68 1.25 
Dingya 22 
表 3 胡麻 产量 形成 模型 产量 构成 因素 法 与 粒 壳 比 法 结果 比较 
Table 3 Comparison of simulation results between yield component factors and grain pod ratio methods 
平均 观测 产量 平均 预测 值 RMSE R 
SA Mean Mean simulated yield (kg:hm “) 
试验 处 理 一 ~ ,= 气 一 . 灿 - 三 所 一 所 LL 三 Pt 
Experiment 。 Treatment measured 产量 构成 因素 法 。 粒 壳 比 法 。 产量 构成 因素 法 粒 壳 比 法 产量 构成 因素 法 “ 粒 壳 比 法 
ye Yield component Grain pod ratio Yield component Grainpod Yield component Grain pod 
(kghm ) factors method method factors method ratio method factors method ratio method 
试验 1 (定西 ) CK 1 059.47 956.13 1 189.50 214.11 391.31 0.82 0.65 
肥料 试验 Yl 1 030.67 760.10 1 242.80 112.14 97.56 0.75 0.62 
Experiment | Y2 1 021.93 867.70 1 204.40 261.35 236.74 0.61 0.76 
(Dingxi): Y3 949.07 861.80 1 107.20 205.45 59.53 0.92 0.66 
站 of RI 938.67 889.00 021.20 71.66 202.61 0.97 0.59 
ertilizer 
R2 914.73 947.90 942.20 24.38 206.76 0.81 0.89 
R3 1 006.80 1 045.70 088.40 27.93 278.39 0.91 0.94 
Fl 1 117.40 1 125.70 132.60 31.48 264.50 0.92 0.66 
F2 1 104.73 1 061.00 240.00 161.11 105.13 0.75 0.64 
F3 1 127.36 1 062.13 274.00 86.03 305.33 0.96 0.58 
试验 中 (定西 ) Tl 1 099.13 814.75 334.88 62.55 20.90 0.87 0.66 
播种 方式 试验 
s T2 934.96 965.00 938.50 213.59 12.62 0.89 0.81 
Experiment || (Dingxi): 
Experiment of T3 1 144.33 1 088.63 295.13 35.68 250.27 0.92 0.71 
sowing method 
试验 川 (定西 ) D1 1 147.67 1 158.33 1 329.17 133.88 269.55 0.61 0.85 
种 植 密度 试验 D2 1 078.72 1 235.50 1 391.00 189.89 268.04 0.63 0.87 
Experiment ll D3 1 240.22 953.50 1 225.33 34.70 398.18 0.81 0.98 
(Dingxi): D4 1 247.83 933.83 1 067.00 95.19 172.83 0.89 0.51 
Experiment of D5 902.89 1 106.17 911.50 203.37 89.65 0.85 0.96 
planting density 
D6 942.33 891.33 959.83 230.54 43.45 0.64 0.58 
D7 931.17 800.17 954.00 274.34 236.92 0.89 0.96 
试验 IV( 榆 中 ) NOPO 1160.33 858.50 1 166.50 48.29 343.72 0.80 0.90 
氮 磷 处 理 试验 NOP1 1 284.33 1 001.00 1 371.50 206.63 21.31 0.87 0.97 
Experiment IV NOP2 998.33 814.50 1 449.50 197.61 62.48 0.88 0.60 
(Yuzhong): Experiment NOP3 953.50 585.00 1 352.00 66.77 307.45 0.66 0.92 
of nitrogen-and N1P0 756.83 541.50 874.50 41.64 86.30 0.61 0.76 
phosphorus level 
NI1Pl1 1 116.67 1 396.50 974.50 146.05 237.40 0.94 0.91 
NI1P2 920.00 971.50 885.50 20.38 358.92 0.90 0.65 
N1P3 946.33 950.50 1 019.50 251.80 38.48 0.95 0.76 
N2P0 833.83 879.50 957.00 173.51 234.73 0.99 0.94 
N2P1 1217.17 1 030.50 1 605.50 208.46 322.86 0.72 0.54 
N2P2 1 669.67 1 796.50 1 195.50 129.92 293.71 0.75 0.94 
N2P3 856.67 648.00 1 422.50 197.34 11.27 0.91 0.77 


CK 为 不 施肥 对 照 ,Y1~Y3 分 别 为 油 汀 施用 量 600 kg:hm“、1 200 kg.hm “和 2 400 kg-hm“, R1~R3 分 别 为 磷酸 二 铵 施用 量 90 kg:hm“、180 kg.hm- 
和 270 kg:hm“, F1~F3 分 别 为 复合 肥 施 用 量 150 kg:hm“、300 kg.hm ”和 450 kg:hm“。T1 为 残 膜 直播 ; T2 为 残 膜 覆 至 春天 ， 播 种 前 揭 残 膜 ， 覆 
盖 新 膜 播种 ; T3 为 播种 前 揭 残 膜 后 直接 播种 .D1~D7 为 7 个 播种 量 处 理 , 分 别 为 3x10° 粒 .hm “、4.5x10° 粒 :hm”、6x10 粒 :hm 7.5x10" 粒 :hm 、 
9x105 粒 .hm ”1.05x107 粒 :hm ”和 1.2x10" 粒 :hm ”。NOP0~N2P3 分 别 代表 3 个 水 平 氮 0 kg(N):hm(N0)、75 kg(N):hm”(N1)、150 kg(N):hm”(N2) 
和 4 个 水 平 磷 0 kg(P;05).hm- (P0)、75 kg(P;05).hm-? (P1) 、150 kg(P20s)-hm“(P2)、225 kg(P205)-hm“(P3) 的 12 种 组 合 。CK represents no 
fertilization as control. Y1, Y2 and Y3 represent sludge fertilizer application rates of 600 kg.hm 一, 1 200 kg.hm and 2 400 kg.hm ,respectively. R1， 
R2 and R3 represent diammonium phosphate fertilizer application rates of 90 kg.hm ,180 kg.hm ”and 270 kg.hm ,respectively. F1, F2 and F3 
represent compound fertilizer application rates of 150 kg:-hm™, 300 kg:hm and 450 kg:hm respectively. T1, T2 and T3 represent sowing methods of 
direct sowing with residual film mulching, sowing with new film mulching after residue film removing in spring, and direct seeding after residual film 
removing in spring, respectively. D1-D7 represent 7 sowing densities of 3x106, 4.5x10°, 6x10°, 7.5x104 9x10°, 1.05x107 and 1.2x107 seeds per hm’. 
NO0OP0 to N2P3 are 12 combinations of 3 levels nitrogen and 4 levels phosphorus. 
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产量 构成 因素 法 
Yield component factors “” 口 
o 粒 壳 比 法 
Grain pod val 


0 500 1000 1500 2000 


© 10 种 肥料 Ten levels of fertilizer 

0 3 种 播种 方式 Three sowing methods 

A 7 个 种 植 密度 Seven planting densities 

x 12 个 氮 磷 水 平 Twelve levels of nitrogen and phosphorus 


模拟 值 Simulated value (kg-hm”) 


0 500 1000 1500 2000 


0 500 1000 1500 2000 
实测 值 
Measured value (kg.hm 7) 

1 胡麻 产量 形成 模型 产量 构成 因素 法 与 粒 壳 比 法 实 
测 值 与 模拟 值 结果 (a)、 产 量 构成 因素 法 (b) 和 粒 壳 比 法 (9 
对 不 同 肥料 、 播 种 方式 、 种 植 密度 和 和 氮 磷 水 平实 测 值 与 模 

拟 值 结果 比较 
Fig. 1 Comparison of simulation results between yield 
component factors and grain pod ratio methods (a), and 
Simulated results by yield component factors method (b) and 
grain pod ratio method (c) with measured values under various 
fertilizers, sowing methods, sowing densities and N and P levels 


3.2.2 ”定西 试验 站 的 检验 比较 
采用 两 种 产量 形成 方法 分 别 对 试验 | ~ 中 定西 


站 的 10 种 肥料 、3 种 播种 方式 和 7 个 种 植 密度 的 产 
量 进行 对 比 检验 。 结 果 表 明 , 产量 构成 因素 法 对 产 
量 的 模拟 效果 优 于 粒 壳 比 法 。 采 用 产量 构成 因素 法 
所 得 不 同 肥料 、 不 同 播种 方式 、 不 同 种 植 密度 产量 
的 模拟 值 RMSE 值 介 于 24.38~274.34 kg-hm“, 平均 
为 133.47 kg.hm-， 而 粒 壳 比 法 的 RMSE 介 于 12.62~ 
398.18 kg:hm“， 平均 为 195.51 kg:hm“。 利用 最 小 二 
乘法 计算 的 决定 系数 及 检验 预测 值 与 观测 值 的 拟 合 
效果 , 产量 构成 因素 法 预测 的 产量 值 与 观测 值 R? 介 
于 0.607 8~0.967 5, 平均 为 0.819 8， 而 粒 壳 比 法 的 
及 介 于 0.514 2~0.978 0, 平均 为 0.743 9。 

3.2.3 ” 榆 中 试验 站 的 检验 对 上 比 

采用 两 种 方法 对 试验 IV 榆 中 试验 站 的 12 个 氮 磷 
处 理 水 平 的 产量 进行 对 上 比 检 验 。 采 用 产量 构成 因素 法 
所 得 产量 模拟 值 的 平均 RMSE 值 为 140.70 kg:hm， 
而 粒 壳 比 法 的 RMSE 平均 为 193.22 kg:hm“。 利用 最 
小 二 乘法 计算 的 决定 系数 R? 验证 预测 值 与 观测 值 的 
拟 合 效果 , 产量 构成 因素 法 预测 的 产量 值 与 观测 值 
及 平均 为 0.832 9, 而 粒 壳 比 法 的 及 平均 为 0.805 8。 
4 讨论 与 结论 
4.1 讨论 

在 产量 构成 三 因素 当中 ,单位 面积 浆果 数 是 胡 
麻生 产 中 对 产量 影响 最 大 的 因素 "变异 最 大 , 不 
同 栽培 条 件 可 相差 1~5 倍 ， 长 期 试验 得 出 基本 上 6 
万 个 戎 果 可 以 获得 0.5 kg 籽粒 。 另 外 两 个 因素 , 每 
果 粒 数 和 粒 重 在 不 同 栽培 条 件 下 ， 对 产量 影响 相对 
较 小 上 变异 幅度 相差 不 超过 1 倍 ， 若 为 同一 品种 ， 
则 一 般 每 果 粒 数 变 化 范围 在 10% 以 内 , 千粒重 在 5% 
以 内 。 而 当 单 位 面积 薄 果 数 达 到 一 定数 量 , 产量 比较 
高 时 ， 则 每 果 粒 数 与 粒 重 对 产量 的 影响 呈 显著 性 中 1。 
因此 , 在 本 模型 中 将 这 3 个 变量 作为 品种 遗传 参数 
针对 不 同 品种 不 同 栽培 方式 , 输入 的 单位 面积 萌 果 
数 、 每 果 粒 数 与 粒 重 都 不 同 , 由 模型 的 验证 结果 得 
出 ， 本 模型 能 较 精 确 地 模拟 胡麻 产量 。 

在 汤 亮 等 中 I 和 张 亚 杰 ? "1 的 研究 中 , 采用 粒 壳 比 
作为 预测 油菜 籽粒 产量 的 方法 , 研究 中 指出 ,油菜 
不 同 品种 的 粒 壳 比 在 不 同 环境 下 有 一 定 变化 ,其 变 
化 受到 库 源 关系 的 限制 ， 机 理 较 为 复杂 ; 在 传 寿 仲 C 
的 研究 中 指出 粒 壳 比 与 作物 产量 呈 极 显著 正 相 关 。 
由 于 影响 胡麻 产量 的 因素 除 遗 传 因素 、 水 分 、 氮 素 、 
播 期 外 ,还 受 其 他 矿物 质 磷 、 钾 等 胁迫 的 影响 ,本 研 
究 中 采用 粒 壳 比 法 模拟 了 胡麻 籽粒 潜在 产量 ， 并 未 
考虑 到 上 述 因 素 以 及 胡麻 粒 壳 上 比 的 变化 规律 对 胡麻 
产量 形成 的 影响 , 因而 ， 导 致 粒 壳 比 法 的 模拟 精度 
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低 于 产量 因素 构成 法 。 
4.2 ”结论 

本 研究 构建 了 基于 品种 遗传 参数 单位 面积 戎 果 
数 、 每 果 粒 数 、 粒 重 与 水 肥 胁 迫 因 子 、 累 积 光 合 速 
率 的 产量 形成 模型 ， 和 基于 粒 壳 比 和 戎 果 干 物质 总 
量 的 产量 形成 模型 。 利 用 定西 试验 站 10 种 肥料 、3 
种 播种 方式 、7 个 种 植 密度 和 榆 中 试验 站 12 个 氮 磷 
处 理 水 平 的 试验 数据 资料 对 两 种 模型 进行 了 较 充分 
的 统计 分 析 与 对 比 检 验 。 结 果 表 明 , 产量 构成 因素 
法 对 产量 的 模拟 效果 优 于 粒 壳 比 法 ,产量 构成 因素 
法 具有 较 高 的 预测 性 和 通用 性 ,在 西北 胡麻 种 植 区 
定西 和 榆 中 地 区 具有 较 好 的 适用 性 。 但 由 于 试验 资 
料 的 限制 ， 本 研究 建立 的 胡麻 产量 形成 模型 并 未 把 
影响 胡麻 产量 的 所 有 因素 都 考虑 进去 ， 如 病虫害 、 
气候 变化 等 ， 以 及 该 模型 在 其 他 胡麻 主 产 区 ， 如 内 
蒙古 、 河 北 、 山 西 等 的 适应 性 验证 都 将 是 今后 进 一 
步 完 善 研究 的 内 容 ， 旨 在 能 探索 出 精确 预测 各 地 区 
胡麻 产量 的 模拟 方法 。 
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